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Motivated by recent experiments observing spin-orbit torque (SOT) acting on the magnetization 
m of a ferromagnetic (F) overlayer on the surface of a three-dimensional topological insulator (TI), 
we investigate the origin of the SOT and the magnetization dynamics in such systems. We predict 
that lateral F/TI bilayers of finite length, sandwiched between two normal metal leads, will generate 
a large antidamping-like SOT per very low charge current injected parallel to the interface. The 
large values of antidamping-like SOT are spatially loealized around the transverse edges of the F 
overlayer. Our analysis is based on adiabatic expansion (to first order in drhidt) of time-dependent 
nonequilibrium Green functions (NEGFs), describing electrons pushed out of equilibrium both by the 
applied bias voltage and by the slow variation of a classical degree of freedom [such as rh{t)]. From it 
we extract formulas for spin torque and charge pumping, which show that they are reciprocal effects 
to each other, as well as Gilbert damping in the presence of SO coupling. The NEGF-based formula 
for SOT naturally splits into four components, determined by their behavior (even or odd) under the 
time and bias voltage reversal. Their complex angular dependence is delineated and employed within 
Landau-Lifshitz-Gilbert simulations of magnetization dynamics in order to demonstrate capability 
of the predicted SOT to efficiently switch m of a perpendicularly magnetized F overlayer. 

PACS numbers: 72.25.Dc, 75.70.Tj, 85.75.-d, 72.10.Bg 


I. INTRODUCTION 


The spin-orbit torque (SOT) is a recently discovered 
phenomenorP® in ferromagnet/heavy-metal (F/HM) 
lateral heterostructures involves unpolarized charge cur¬ 
rent injected parallel to the F/HM interface induces 
switching or steady-state precessioiP of magnetization 
in the F overlayer. Unlike conventional spin-transfer 
torque (STT) in spin valves and magnetic tunnel junc¬ 
tion (MTJs),^^^ where one F layer acts as spin-polarizer 
of electrons that transfer torque to the second F layer 
when its free magnetization is noncollinear to the direc¬ 
tion of incoming spins, heterostructures exhibiting SOT 
use a single F layer. Thus, in F/HM bilayers, spin-orbit 
coupling (SOC) at the interface or in the bulk of the 
HM layer is crucial to spin -pola rized injected current 
via th e Ede lstein effect (EE)I®^ or the spin Hall effect 
(SHE),f^^^i^ respectively. 

The SOT offers potentially more efficient magnetiza¬ 
tion switching than achieved by using MTJs underlying 
present STT-magnetic random access memories (STT- 
MRAM).C^ Thus, substantial experimental and theo¬ 
retical efforts have been focused on identifying physi¬ 
cal mechanisms behind SOT whose understanding would 
pave the way to maximize its value by using optimal 
materials combinations. For example, very recent ex¬ 
periment Eave replaced HM with three-dimensional 
topological insulators (3D TIs).II3The TIs enhanc^I^J^ 
(by a factor hvpjaR, where vp is the Fermi velocit y on 
the surface of TI and ap is the Rashba SOC strengtfP^ 
at the F/HM interface) the transverse nonequilibrium 
spin density driven by the longitudinal charge current, 
whi ch is responsible for the Idiige ^Id-like SOT compo- 
nentpSESl observed experiment ally. 


FIG. 1. (Golor online) Schematic view of F/TI lateral bilayer 
operated by SOT. The F over layer has finite length and rh 
is the unit vector along its free magnetization. The TI layer 
is attached to two N leads which are semi-infinite in the x- 
direction and terminate into macroscopic reservoirs. We also 
assume that F and TI layers, as well as N leads, are infinite in 
the y-direction. The unpolarized charge current is injected by 
the electrochemical potential difference between the left and 
the right macroscopic reservoirs which sets the bias voltage, 
Ml — I^R = eVb . We mention that the results do not change if 
the TI surface is covered by the F overlayer partially or fully. 


Furthermore, recent experiments have also observed 
antidamping-like SOT in F/TI heterostructures with 
surprisingly large figure of merit (i.e., antidamping 
torqu^er unit applied charge current density) that sur- 
passes^^^^ those measured in a variety of F/HM het¬ 
erostructures. This component competes against the 
Gilbert damping which tries to restore magnetization 
to equilibrium, and its large figure of merit is, there¬ 
fore, of particular importance for increasing efficiency of 
magnetization switching. Theoretical understanding of 
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the physical origin of antidamping-like SOT is crucial to 
resolve the key challenge for anticipated applications of 
SOT generated by TIs—demonstration of magnetization 
switching of the F overlayer at room temperature (thus 
far, magnetization switc hing has been demonstrated only 
at cryogenic temperature^^ . 

However, the microscopic mechanism behind its large 
magnitude^HH and ability to efficiently (i.e., using as lit¬ 
tle dc current density as possible) switch magnetizat ioiP^ 
remains under scrutiny. For example, TI samples used 
in these experiments are often unintentionally doped, so 
that bulk charge carriers can generate antidamping-like 
SOT via rather lar gS HE (but not sufficient to explain 
all reported value^^^^^ . The simplistic picture,!^ in 
which electrons spin-polarized by the EE diffuse into the 
E overlayei^i^ to deposit spin angular momentum within 
it, cannot operate in technologically relevant E overlayers 
of ^ 1 nm th icknes or explain complex angular depen- 
denc^^E^EH typic ally o bserved for SOT. The Berry cur¬ 
vature mechanisnPSI^ for antidamping-like SOT applied 
to lateral E/TI heterostructures predicts its peculiar de¬ 
pendence on the magnetization orientation,!^ vanishing 
when magnetization fa is parallel to the E/TI interface. 
This feature has thus far not been observed experimen- 
tally,II^and, furthermore, it makes such antidamping-like 
SOT less efhcienfP^ (by requiring larger injected currents 
to initiate magnetization switching) than standard SHE- 
driverP^ antidamping-like SOT. 

We note that the recent experimentaP^^ll theo- 
reticaP^^ studies of SOT in lateral E/TI bilayer have 
focused on the geometry where an infinite E ov erlaye r 
covers an infinite TI layer. Moreover, they assume^^l^ 
purely two-dimensional transport where only the top sur¬ 
face of the TI layer is explicitly taken into account by the 
low-energy effective (Dirac) Hamiltonian supplemented 
by the Zeeman term due to the magnetic proximity ef¬ 
fect. On the other hand, transport in realistic Tl-based 
heterostructures is always three-dimensional, with unpo¬ 
larized electrons being injected from normal metal con¬ 
tacts, reflected from the E/TI edge to flow along the sur¬ 
face of the TI in the ^z-plane and then along the bottom 
TI surface in Eig. In fact, electrons also flow within a 
thin layer (of thickness < 2 nm in Bi 2 Se 3 as the prototyp¬ 
ical TI material) underneath the top and bottom surfaces 
due to top and bottom metallic surfaces of the TI doping 
the bulk via evanescent wave functions.^ Therefore, in 
this stu^ we consider more realistic and experimentally 
relevanfP^ E/TI bilayer geometries, illustrated in Eig. [h 
where the TI layer of finite length and finite thick¬ 
ness is (partially or fully) covered by the E overlayer 
of length . The two semi-infinite ideal N leads are di¬ 
rectly attached to the TI layer, we should mention that 
the result does not depend on the length of TI layer that 
is covered by the EM. 

Our principal results are twofold and are summarized 
as follows: 

(z) Theoretical prediction for SOT: We predict that the 
geometry in Eig. will generate large antidamping-like 


SOT per low injected charge current. By studying spatial 
dependence of the SOT (see Eig. |^, we show that in a 
clean EM/TI interface the electrons exert anti-damping 
torque on the EM as they enter into the interface and un¬ 
less interfacial roughness or impurities are included the 
torque remains mainly concentrated around the edge of 
the interface. Although the exact results show strong 
nonperturbative features, based on second order pertur¬ 
bation we present two different interpretations showing 
that the origin of the antidamping SOT relies on the 
spin-flip reflection of the chiral electrons injected into 
the EM/TI interface. Its strong angular dependence (see 
Fig.§, i.e., dependence on the magnetization direction 
m, offers a unique signature that can be used to distin¬ 
guish it from other possible physical mechanisms. By 
numerically solving the Landau-Lifshitz-Gilbert (LEG) 
equation in the macrospin approximation, we demon¬ 
strate (see Eigs.j^and that the obtained SOT is ca¬ 
pable of switching of a single domain magnetization of a 
perpendicularly magnetized E overlayer with bias voltage 
in the oder of the Magneto-Crystaline Anisotropy (MCA) 
energy. 

(a) Theoretical formalism for SOT: The w idely used 
quantum (such as the Kubo formulaP ^^^ I ^^ !) and semi- 
classical (such as the Boltzmann equatioiP^ transport 
approaches to SOT are tailored for geometries where an 
infinite E layer covers an infinite TI or HM layer. Due to 
translational invariance, the nonequilibrium spin density 
S induced by the EE on the surface of TI or HM layer 
has uniform orientation S = (0,5'y,0) [in the coordinate 
system in Eig. [^, which then provides reference direc¬ 
tion for defining field-like, Tfm x and antidamping- 
like, Tadfi^ X (m X ^), components of SOT. In order to 
analyze spatial dependence of SOT in the device geome¬ 
try of Eig.[^ while not assuming anything a priori about 
the orientation of field-like and antidamping-like compo¬ 
nents of SOT, we employ adiabatic expansiorP^ of tinie- 
dependent nonequilibrium Green functions (NEGEsp^^ 
to derive formulas for torque, charge pumpin^^^® and 
Gilbert dampin^^ in the presence of SOC. The NEGE- 
based formula for SOT naturally splits into four compo¬ 
nents, determined by their behavior (even or odd) under 
the time and bias voltage reversal. This gives us a general 
framework in quantum mechanics to analyze the dissi¬ 
pative (antidamping-like) and nondissipative (field-like) 
force (torque) vector fields for a set of canonical variables 
(magnetization directions). Their angular (see Eig. 
and spatial (see Eig. dependence shows that although 
field-like and antidamping-like SOTs are predominantly 
along the fhxy and fax {fax y) directions, respectively, 
they are not uniform and can exhibit significant devia¬ 
tion from the trivial angular dependence defined by these 
cross products [see Eig. ih)]. 

The paper is organized as follows. In Sec.|llj we present 
the adiabatic expansion of time-dependent NEGEs, in a 
representation that is alternative to Wigner representa- 
tiojjm (usually employed for this type of derivatiorP^ , 
and derive expressions for torque, charge pumping and 
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Gilbert damping. In Sec. |ml we decompose the NEGF- 
based expression for SOT into four components, deter¬ 
mined by their behavior (even or odd) under the time and 
bias voltage reversal, and investigate their angular de¬ 
pendence. Section [TV| discusses the angular dependence 
of the zero-bias transmission function which identifies the 
magnetization directions at which substantial reflection 
occurs. In Sec. [V| we study spatial dependence of SOT 
components and discuss their physical origin. Section [VT| 
presents LEG simulations of magnetization dynamics in 
the presence of predicted SOT, as well as a switching 
phase diagram of the magnetization state as a function 
of the in-plane external magnetic field and SOT. We con¬ 
clude in Sec. EUl 


II. THEORETICAL FORMALISM 


We first describe the time-dependent Hamiltonian 
model, H(t) = Ho+U(t), of the lateral F/TI heterostruc¬ 
ture in Fig. B Here Hq is the minimal tight-binding 
model for 3D TIs like Bi 2 Se 3 on a cubic lattice of spacing 
a with four orbitals per site.^^ The thickness, = 8a 
of the TI layer is sufficient to prevent hybridization be¬ 
tween its top and bottom metallic surface statcs.^^^ The 
time-dependent potential 

U(t) = -Asurflm^(^) • d'/2, (1) 


depends on time through the magnetization of the F over¬ 
layer which acts as the slowly varying classical degree of 
freedom. Here m{t) is the unit vector along the direc¬ 
tion of magnetization, Agurf = 0.28 eV is the proximity 
induced exchange-field term and is a diagonal matrix 
with elements equal to unity for sites within the F/TI 
contact region in Fig. and zero elsewhere. The semi¬ 
infinite ideal N leads i n Fig . are taken into account 
through the self-energieP^®^ computed for a tight- 
binding model with one spin-degenerate orbital per site. 
The details of how to properly couple to Hq, while 

taking into account that the spin operators for electrons 
on the Bi and Se sublattices of the TI are inequivalent ,1^ 
can be found in Ref. sni 

Within the NEGF formalisnP^^^ the advanced 
and lesser GFs matrix elements of the tight-binding 
Hamiltonian, Hq, are defined by = 

-i@{t - and G5, t') = 

^(4o's'(^0cios(i)), respectively. Here, cJos ^he 

creation (annihilation) operator for an electron on site, i, 
with orbital, o, and spin 5, respectively, (...) denotes the 
nonequilibrium statistical average, and h = 1 to simplify 
the notation. These GFs are the matrix elements of the 
corresponding matrices G and used throughout the 
text. 

Under stationary conditions, the two GFs depend on 
the difference of the time arguments, t — t', and can 
be Fourier transformed to energy. In the strictly adi¬ 
abatic limit one can emplo}!^ the same retarded GF, 


Gt(F) = [E — H(t) — Hl — as under stationary 

conditions, but where the GF depends parametrically on 
time (denoted by the subscript t) and is computed for the 
frozen-in-time configuration of U(t). However, even for 
slow evolution of m{t) corrections^ to the adiabatic GF 
are needed to describe dissipation effects such as Gilbert 
damping or the charge current which can be pumped^ 
by the dynamics of m{t). 

The so-called adiabatic expansion, which yields correc¬ 
tions beyond the strictly adiabatic limit, is traditionally 
performed using the Wigner representatiorP^ in which 
the fast and slow time scales are easily identifiable.l^The 
slow motion implies that the NEGFs vary slowly with the 
central time tc = (t + t')/2 while they change fast with 
the relative time U = t — t' . By expanding the Wigner 
transformation of NEGFs 

^c) = J dtr ~ 5 

( 2 ) 

in the central time tc while keeping only terms contain¬ 
ing first-order derivatives d/dtc (due to the slow varia¬ 
tion with tc) gives the first-order correction beyond the 
strictly adiabatic limit This route requires to han¬ 
dle complicated expressions resulting from the Wigner 
transform applied to convolutions of the type G(ti, ^ 2 ) = 

J dts Gi(ti, t3)G2(t3, ^ 2 ). 

Here we provide an alternative derivation of the first- 
order nonadiabatic correction. Namely, we consider t 
(observation time) and t — t' (relative time) as the nat¬ 
ural variables to describe the time evolution of NEGFs 
and then perform the following Fourier transfornP^ 

G{t,t') = p (3) 


The standard equations of motion for G{t ,t') and 
G^(t,t') are cumbersome to manipulat^^^® or solve 
numerically,!^ so they are usually transformed to some 
other representation.!^ Here we replace G(t,t') in the 
standard equations of motion with the rhs of Eq. (§ to 
arrive at: 




G{E,t) = l, 

(4) 


and 

G<{t,t') = J ^G{E,t)'S<{E)G\E,t')e^^^^-^'\ (5) 

For the two-terminal h eteros tructure in Fig. in the 
elastic transport regime,!^^!^ ^(^) = ^l{E) + T,ji{E) 
and J:<{E) = ifL{E)TL{E) + ifR{E)TR{E), where 
El,r = ~ Fermi-Dirac distribution 

functions of electrons in the macroscopic reservoirs into 
which the left and right N leads terminate are ,r{E) = 
f{E — iil.r)^ where the difference between the electro¬ 
chemical potentials, iil.r = Ep -h cVl^r, defines the bias 
voltage eVb = I^l - I^r^ 
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Using the following identity 


^ iGr„Gt - (G - Gt) + i A (^G + G , 


a=L,R 


(6) 


the lesser GF to first order in small d\J{t)/dt and in the low bias Vb regime (i.e., the linear-response transport regime) 
can be expressed as, 


j ^ l[G{E,t)-G\E,t)]fiE)+ feVaGtT^Gl+ifGt 

^ \ a=L,R 


dt 


( 7 ) 


where, the first term corresponds to the density matrix of the equilibrium electrons occupying the time dependent 
single particle states, while the second and third terms describe the density matrix of the excited (nonequilibrium) 
electrons occupying the states close to the Fermi energy due to the bias voltage and time dependent term in the 
Hamiltonian, respectively. The retarded GF in the first term in Eq. 0 expanded to first order in 9U(t) /dt is of the 
form 


G{E,t)^Gt + i^^^Gt. 

The lesser GF determines the time-dependent nonequilibrium density matrix 


(8) 


pit) = -G^{t,t), 


( 9 ) 


from which we determine the time-dependent expectation values of physical observables, A{t) = Tr[p(t)A]. In 
particular, the relevant quantities for the heterostructure in Fig. are the charge current 


/ dW 

— IV {r„(E)G< (i, t) + UiE)T^{E) [G{E, t) - Gt {E, t)]} 


27r 


IdEf{E) I J2iyp - Ve.)T^^iE) 


Asurf drUj 


2e 


dt 


( 10 ) 


and the spin density 




dE 

27r 


Y,UE)R%E) -^Y.f{E)^Tt{E) 


( 11 ) 


where a, 13 G {L,R} and i,j G {x,y,z}. 


The “trace-formulas” in Eqs. ( |10[ ) and ( |ll ) 
= Tr 


= Tr 
^“(E) = Tr 


r^GtT^jGl 

^mf^iGlT aGt 
imC’'iGtTaGl 


( 12 a) 

( 12 b) 

( 12 c) 


T^^{E) = Tr \lmCT^{Gl - G,)l^t 7 ,(G, - Gl)l(,12d) 


determine charge currenfP^ due to V 5 , charge current 
pumpecP^ by the dynamics of m{t) in the presence of 
SOG, the spin torque, and Gilbert damping tensor, re¬ 
spectively. In the expression for the spin density we ig¬ 
nore the antisymmetric part of which corresponds to 
the renormalization of the precession frequency of mag¬ 


netization dynamics.l^Note that in Eq. (p!^(b) and 
its time-reversal in Eq. (12)(c) reveal a reciprocal 
relation between charge pumping by magnetization dy¬ 
namics in the presence of SOG and current-driven SOT 
at each instant of time t. 

The spin density in Eq. 0 enters into the LL equa¬ 
tion for the magnetization dynamics 


dm 

dt 


^surf 


m X s{t). 


(13) 


The first term in Eq. 0 generates the spin torque in 
Eq. (13) which has three contributions: (z) an equilib¬ 
rium component responsible for the interlayer exchange 
interaction in the presence of a second E layer and/or 
magneto-crystalline anisotropy (MGA) in the presence of 
SOG; (a) a bias-induced field-like torque modifying the 
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equilibrium interlayer exchange and MCA fields; and {in) 
a damping (antidamping)-like torque describing angular 
momentum loss (gain) due to the flux of Fermi surface 
electrons. 


III. ANGULAR DEPENDENCE OF SOT 
COMPONENTS 

In order to understand the different contributions to 
T*^(F^), we decompose it into even (e) or odd (o) terms 
under time-reversal Gt gJ : 

Ti^{E) = [T^{E) + T^\E)]/2, (14) 

Ti^{E) = [T^{E) - T^\E)]/2. (15) 

Since '^o^Tl^{E) = 0, the contribution of the odd com¬ 
ponent to the equilibrium spin density in Eq. van¬ 
ishes identically, while its nonzero values appear only for 
E within the bias window around Ep. This motivates 
further splitting of T'^^{E) into four components for the 
case of two-terminal devices 


, Ti\E) + F^{E) 

= 2 

(16a) 

, Ti^{E)-Ti^{E) 

= 2 

(16b) 

where the first and second subscripts denote their be¬ 
havior (even or odd) under bias reversal V 5 —Vh and 

time reversal, respectively. The corresponding four com¬ 
ponents of torque are determined by 

fe,. = 1 ^ [fdE) + fn{E)%AE), 

(16c) 

fo,u = ifUE) - ME)]toAE), 

(16d) 

where the energy-resolved torque is given by 


F,{E) = -AmxFFE), 

(17) 


and E {e, o}. The terms To^o and To^e are non-zero 
only in nonequilibrium driven by 14 7 ^ 0, and depend on 
electronic states in the bias voltage window around the 
Fermi energy (or on the Fermi surface states in the linear- 
response regime where integrals are avoided by multiply¬ 
ing integrand by eV^). The term Te,o = 0 is zero, while 
Te,e is nonzero also in equilibrium and, therefore, depends 
on all occupied electronic states. 

Figures l^a-c) show the net vector field (summed over 
all sites of the F overlayer) of T^^y{EF) at zero temper¬ 
ature for different directions of m on the unit sphere. 
Their angular behavior reveals that: (z) Tq,© shown in 
Fig. [^a) is the field-like SOT generated by the EE, with 
predominant orientation along the m x ^-direction; (ii) 
ro,e in Fig. |^b) is the antidamping-like SOT with pre¬ 
dominant orientation along the m x {m x ^)-direction; 
and (Hi) re,e ia Eig. |^c) along the rh x z direction is 


the field-like component whose angular dependence be¬ 
haves approximately as 2(m • z)|m x = sin(26>) typical 
for torque components generated by the MCA field.^ 


The corresponding angular dependence of the net 
Tq o(^ G {x^y^z}), |Fo^, and rj g along the solid trajecto¬ 
ries shown in Eigs. [^a-c) are plotted in Eigs. id-f), 
respectively. We find that the magnitude of the max¬ 
imum antidamping-like SOT is about a factor of four 
larger than that of the field-like SOT. Additionally, the 
field-like SOT peaks when the magnetization is in-plane. 
In contrast, the antidamping-like SOT peaks when the 
magnetization is out of plane, which can be attributed to 
the gap opening of the Dirac cone which in turn enhances 
the reflection (see Sec. IV) at the lateral boundaries of 
the E over layer. 


Note that the magnitude of the net SOT components 
shown in Eigs. is- i) exhibits strong angular dependence 
because of the large SOC on the TI surface similar to that 
found in E/HM heterostructures wh en the Rashba SOC 
at the interface is sufficiently strong.1^^1^ In particular, 
the significant deviation of the angular dependence of the 
antidamping-like SOT from the trivial |m x (m x ^)| be¬ 
havior in the limit of m ^ y, shown quantitatively in 
Fig. Hh), indicates its nonperturbative variation with 
respect to the magnetization direction. A similar nonper¬ 
turbative angular behavior (i.e., strong deviation from 
the standard oc sin^ 0 dependence on the precession cone 
angle 0) has been found for adiabatic charge pumping 
fro m a p recessing E overlayer attached to the edge of 2D 
rpjj47l48] surface of 3D TIs.^ 


IV. ANGULAR DEPENDENCE OF 
TRANSMISSION FUNCTION 


Figure]^ shows the transmission function T^^{Ef) in 
Eq. ( |12a[ ) for the heterostructure in Eig. [^versus the ori¬ 
entation of rh on the unit sphere. We find that the charge 
current determined by T^^{E) is smallesfP^ when m \\ z 
or m II X. This is due to the reflection of Dirac elec¬ 
trons on the TI surface from the lateral boundaries of 
the E overlayer. Underneath the E overlayer, the ex¬ 
change field, — Asurf^ • d^/ 2 , induced by the magnetic 
proximity effect is superimposed on the Dirac cone sur¬ 
face dispersion. This opens an energy gap Agurf when 
m II z (or smaller gap Ag^rf cos 6 > for ruz 7 ^ 0) at the DP 
of the TI region underneath the E overlayer, which in 
turn gives rise to strong electronic reflection when rh \\ z. 
Eor m II X, there is no energy gap at the DP and the 
Dirac cone effectively shifts away from the center of the 
Brillouin zone due to proximity exchange field. Neverthe¬ 
less electrons polarized by the EE along the ^-axis reflect 
from the magnetization pointing along the x-axis. 
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(a) 





FIG. 2. (Color online) (a)-(c) The vector field of SOT components defined by Eq. (16), for different directions of 

m on the unit sphere. The angular behavior in (a) and (c) shows that To,o and Te,e behave as held-like torques, while that in 
(b) shows that To,e behaves as antidamping-like torque. Cartesian components [(d)-(f)] and magnitude of [(g)“(i)] along 
the trajectories denoted by solid lines in the corresponding panels (a)-(c). The magnitude of is divided by: |m x y| in (g); 
\m X (m X y)\ in (h) and 2(m • z)\m x zj in (i). 


V. SPATIAL DEPENDENCE OF SOT 
COMPONENTS AND PHYSICAL ORIGIN OF 
ANTIDAMPING-LIKE SOT 


Figure |^a) demonstrates that large values of 
antidamping-like SOT from Fig. are spatially localized 
around the transverse edges of the F overlayer, for Fermi 
energy inside and outside of the surface state gap induced 
by out-of-plane magnetic exchange coupling. While con¬ 
ductance in this system is close to zero at the Dirac point, 
we observe that the anti-damping torque does not de¬ 
pend strongly on the Fermi energy. This suggests a high 
efficiency of SOT per injected current for the Fermi en¬ 
ergy close to the Dirac point. For the field-like SOT 
in Fig. Sb) we see the torque independent of the co¬ 
ordinate in the entire F/TI contact region, as expected 
from the phenomenology of the EE. In Fig. [^c) we plot 
the contribution of the Fermi energy electrons to the 
FM/TI interface induced MCA field. Even though to¬ 
tal = 0, its spatially-resolved value plotted in 

Fig. md) is nonzero which can be removed by perform¬ 
ing a proper gauge transformation. 

To understand the origin of the anti-damping SOT let 
us find an expression for the average of SOT around a 


fixed axes m = mo-\-{rh± cos(0)+mo xm_L sin(0))(5^ with 
small cone angle SO. By applying a rotation operator we 
can align the fixed axes along 2 :-axis such that mo = Cz 
and m± = Cx- In this case we have 

(r-), = j I f f ( 18 ) 

a 

where, and ^ refers to the imaginary 

part. Expanding the GFs in Eq.pl^b) to the lowest order 
with respect to S0e~'^^^ we obtain. 





l^Cr ( 19 ) 


Plugging this expression into Eq.( 18), and using the iden¬ 
tity, Gt — gJ = i in linear bias voltage 

regime we obtain. 


( r ^)0 = 




surf*^^ rp r tt 1 






16 


where, ~ corresponds to the 

density matrix inside the F overlayer for the electrons 
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FIG. 3. (C olor online) The transmission function 
in Eq. (12a) for two-terminal heterostructure shown in Fig.[^ 
at different directions of magnetization ifi on the unit sphere. 
The Fermi energy is set at Ff = 3.1 eV (which is 0.1 eV above 
the DP). 



Longitudinal Coordinate (a) 

FIG. 4. (Color online) (a)-(d) Spatial dependence of SOT 
components, t^ j^(Ff) (i G {x,y,z} and /x, G {e, o}), per 
unit length, for Fermi energy inside the magnetization in¬ 
duced gap around Dirac point {Ef = S.leV) and outside 
the gap {Ef = 3.3eV), for m \\ z in Fig. Their physical 
meaning is explained in Fig. The range of x-coordinate 
corresponds to the length = 40a of the F overlayer, while 
the results are independent of the length of the TI layer un¬ 
derneath. 


(holes) at the Fermi surface being injected from the lead 
a {P ^ a). The electron-hole analogy can be understood 
by defining the hole density matrix, iG ^, from the idem 
tity -i{G< -G>)= 2?j(G) = Pa + PP- % con- 
sidering left-lead induced holes instead of right-lead in¬ 
duced electrons, we can interpret Eq.(2^ as spin-resolved 
electron-hole recombination rate, where opposite spins 
have opposite contributions to the antidamping-like 
SOT. This picture focuses on the energy anti-dissipative 
aspect of the phenomena and, since P^{p^r) cor- 


eV,=0 eV,=-4.5 A^ca 

FIG. 5. (Color online) SOT-induced magnetization trajecto¬ 
ries m{t) under different 14 and = 0. Higher color inten¬ 
sity denotes denser bundle of trajectories which start from all 
possible initial conditions rh{t = 0) on the unit sphere. Solid 
curves show examples of magnetization trajectories, while the 
white circles denote attractors of trajectories. 




responds to the spin-cr right (left) moving electrons, 
Eq. (20) suggests that spin-momentum locking natu¬ 
rally has a significant effect on the enhancement of 
the antidamping-like SOT magnitude. In particular, 
in the case of F/TI interface, the enhancement of the 
antidamping-like SOT occurs when the spin-up/down is 
along the ^-axis (mo || y) which is the spin-polarization 
direction of electrons passing through the surface of the 
TI induced by the EE. Additionally, in this case the 
antidamping-like SOT gets smaller away from the F/TI 
transverse edge because the contribution of both of the 
leads to the spin density become identical. Therefore the 
anti-damping torque in this case is more localized around 
the edge. This effect is more significant when the magne¬ 
tization is out of the plane and the Fermi energy is inside 
the Asurf cos 0 gap on the TI surface. 

A alternative interpretation of the results can be 
achieved by considering Gt — gJ = Gjrc^Gt. In 
this case, the average of the antidamping-like SOT is ex¬ 
pressed by 


{T% = ^Tr[Tli-Tli], 


( 21 ) 


where the F overlayer induced spin-flip transmission ma¬ 
trix is defined as 


and 


rrit-l _ A-/-74 

“^a/3 ~ V^a/3/ ^a/3’ 


^al3 ~ 9 




( 22 ) 


(23) 


Although Eq. (21) is obtained from perturbative con¬ 


siderations, it looks identical to the Eq. (8) of Ref. [48] 
where a spin-flip reflection mechanism at the edge of the 
F/2D-TI interface was recognized to be responsible for 
the giant charge pumping (i.e., anti-damping torque) ob¬ 
served in the numerical simulation.SSlEq. ( [2^ describes a 
transmission event in which electrons injected from lead 
(a, get spin-flipped (from up to down) by the FM and 
then transmit to the lead p. The path of the electrons 
describing this process is shown in Fig. From the k- 
resolved results of the anti-damping torque (not shown 






















here) we observe that while for the in-plane magnetiza¬ 
tion electrons moving in the same transverse direction 
(same sign for ky) on both left and right edges of the 
FM/TI interface contribute to the torque, in the case of 
out-of-plane magnetization for the left (right) edge of the 
interface the local anti-damping torque is induced mostly 
by the electrons with ky >0 {ky < 0 ). 

It is worth mentioning that due to nonperturbative na¬ 
ture of the SOT induced by the chiral electrons, the ap¬ 
proximation presented in this section which can as well 
be obtained from the self energy corresponding to the 
vacuum polarization Feynman diagrams of the electron- 
magnon coupled systeni^, does not capture the phenom¬ 
ena accurately. This is evident in the angular dependence 
of the anti-damping torque which in the current section 
is considered up to second order effect {S0‘^), while the 
divergence-like behavior in Figsj^h) suggest a linear de¬ 
pendence when the magnetization direction is close to the 
^-axis. This signifies the importance of the higher order 
terms with respect to 60 that can not be ignored. The ap¬ 
proximation presented in this section also suggests that 
blocking the lower surface leads to the reduction of the 
anti-damping torque. However, in this case an electron 
experiences multiple spin-flip reflections before transmit¬ 
ting to the next lead and in fact it turns out that the ex¬ 
act results stay intact even if the lower surface is blocked. 
This is similar to the conclusion made in Ref. [48] which 
shows the redundancy of blocking the lower edge of the 
2D-TI to obtain a nonzero pumped charge current from 
processing FM as proposed in Ref. iTTl 

Although spin-momentum locking of the surface state 
of the TI resembles the 2D Rashba plane, in the case of 
TI surface state the cones with opposite spin-momentum 
locking reside on opposite surface sides of the TI slab 
while in the case of a Rashba plane they are only sepa¬ 
rated by the SOC energy. This means one can expect a 
smaller SOT for a FM on top of a 2 D Rashba plane due 
to cancellation of the effects of the two circles with op¬ 
posite spin-momentum locking, where the nonzero anti¬ 
damping torque originates from the electron-hole asym¬ 
metry. 


VI. LLG SIMULATIONS OF MAGNETIZATION 
DYNAMICS IN THE PRESENCE OF SOT 


In order to investigate ability of predicted 
antidamping-like SOT to switch the magnetization 
direction of a perpendicularly magnetized F overlayer in 
the geometry of Fig. we study magnetization dynamics 
in the macrospin approximation by numerically solving 
LLG equation at zero temperature supplemented by 
SOT components analyzed in Sec. in 


dm 

dt 


= T[^o,e(w, Ef) + To,o{m, EF)]eVb + x m 

ZTT 


+m X 


. dm 

c(m) ■ 


+ (m • z){m X i)AMCA- (24) 


6 

■R 




FIG. 6. (Color online) Phase diagram of the magnetization 
state in lateral F/TI heterostructure from Fig.[^as a function 
of an in-plane external magnetic || x and Vb (i.e., SOT 
oc Vb). Thick arrows on each of the panels (a)-(d) show the 
direction of sweeping of or Vb parameter. The small¬ 
ness of central hysteretic region along the U-axis, enclosed 
by white dashed line in panel (b) and (d), shows that low 
currents are required to switch magnetization from > 0 
to mz < 0 stable states. 


Here 7 is the gyromagnetic ratio, cx{m)ij = 
is the dimensionless Gilbert damp¬ 
ing tensor, and Amca = ^mca ^ \fe,e\/\{rn ’ z){m x z)\, 
where represents the intrinsic MCA energy of the 

FM. We solve Eq. (24) by assuming that the Gilbert 
damping is a constantjits dependence on m is relegated 
to future studies) and ignore the dependence of Amca 
on m and V 5 while retaining its out-of-plane direction. 

Figure shows the magnetization trajectories for all 
possible initial conditions m{t = 0 ) on the unit sphere 
under different V^. At W = 0 , the two attractors are 
located as the north and south poles of the sphere. At 
finite Vb, the attractors shift away from the poles along 
the z-axis within the xz-plane, while additional attractor 
appears on the positive (negative) ^-axis under negative 
(positive) Vb. Note that the applied bias voltage Vb drives 
dc current and SOT proportional to it in the assumed 
linear-response transport regime. _ 

Figure shows the commonly constructecP^^^^^^^ 
phase diagram of the magnetization state in the pres¬ 
ence of an external in-plane magnetic field || x and 
the applied bias voltage Vb (i.e., SOT (x Vb). The thick 


denote the direction of the 
a) §b)] we increase [de¬ 


arrows in each panel of Fig. 6 
sweeping variable—in Fig. 6 
crease] Vb slowly in time, and similarly in Fig. [^c) j^d)] 
we increase [decrease] the external magnetic field gradu¬ 
ally. The size of hysteretic region in the center of these 
diagrams, enclosed by white dashed line in FigsJ^b) and 
§d), measures the efficiency of switchingSince 
this region, where both magnetization states mz > 0 and 
mz < 0 are allowed, is relatively small in Figs. [^a) and 
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ib), magnetization can be switched by low and 

small Vb (or, equivalently, small injected dc current), akin 
to the phase diagrams observed in recent experiments.^ 
Although we considered the FM a single domain, the 
fact that the anti-damping component of the SOT is 
mainly peaked around the edge of the FM/TI interface 
suggests that it is be more feasible in realistic cases to 
have the local magnetic moments at the edge of the FM 
switch first and then the total magnetization switches 
by the propagation of t he do main walls formed at the 
edge throughout the Therefore, a micromag- 

netic simulation of the system is required to investigate 
switching phenomena in large size systems which we rel¬ 
egate to future works. 


VII. CONCLUSIONS 

In conclusion, by pe rform ing adiabatic expansion of 
time-dependent NEGFs,l^^^^we have developed a frame¬ 
work which yields formulas for spin torque and charge 
pumping as reciprocal effects to each other connected by 
time-reversal, as well as Gilbert damping due to SOG. It 
also introduces a novel way to separate the SOT com¬ 
ponents, based on their behavior (even or odd) under 
time and bias voltage reversal, and can be applied to 
arbitrary systems dealing with classical degrees of free¬ 
dom cou ple d to electrons out of equilibrium. For the 
geometry^ proposed in Fig. where the F overlayer 


covers (either partially or fully) the top surface of the TI 
layer, we predict that low charge current flowing solely 
on the surface of TI will induce antidamping-like SOT on 
the F overlayer via a physical mechanism that requires 
spin-momentum locking on the surface of TIs—spin-flip 
reflection at the lateral edges of a ferromagnetic island 
introduced by magnetic proximity effect onto the TI sur¬ 
face. This mechanism has been overlooked in efforts to 
understand why SO-coupled interface alone (i.e., in the 
absence of SHE current from the bulk of SO-coupled non¬ 
ferromagnetic materials) can generate antidamping-like 
SOT, where other explored mechanisms have included 
spin-dependent impurity sc atteri ng at the interface,!^ 
Berry curvature mechanism,as well as their com¬ 
bination.!^ 

The key feature for connecting experimentally ob¬ 
served SOT and other related phenomena in E /TI het- 
erostructures (such as spin-to-charge conversiorP^^^^l^ 
to theoretical predictions is their dependenc^^^^ on the 
magnetization direction. The antidamping-like SOT pre¬ 
dicted in our study exhibits complex angular dependence, 
exhibiting “nonperturbative” change with the magnetiza¬ 
tion direction in Eig.|^h), which should make it possible 
to easily differentiate it from other competing physical 
mechanisms. 
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